Many-body Electronic Structure of Metallic ex-Uranium 
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We present results for the electronic structure of a uranium using a recently developed quasipar- 
ticle self-consistent GW method (QSGW). This is the first time that the /-orbital electron-electron 
interactions in an actinide has been treated by a first-principles method beyond the level of the 
generalized gradient approximation (GGA) to the local density approximation (LDA). We show 
that the QSGW 7 approximation predicts an /-level shift upwards of about 0.5 eV with respect to 
the other metallic s-d states and that there is a significant /-band narrowing when compared to 
LDA band-structure results. Nonetheless, because of the overall low /-electron occupation number 
in uranium, ground-state properties and the occupied band structure around the Fermi energy is 
not significantly affected. The correlations predominate in the unoccupied part of the / states. This 
provides the first formal justification for the success of LDA and GGA calculations in describing the 
ground-state properties of this material. 

PACS numbers: 71.15-m,71.10-w,71.20-Eh 
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It has long been recognized that 5/ electron-electron 
correlations play an important role in the light actinidcs 
0, Q, becoming increasingly significant as one moves 
across this series and the atomic number Z increases. 
This culminates in Pu, which has many extreme phys- 
ical properties that are driven by these correlations [3|, 
such as the large volume expansion for the a to 5 phase 
transformation 0, HJ . What is less clear is the role of cor- 
relations for Z's less than Pu. Uranium stands at a kind 
of threshold in this regard. Experimentally, the pure ma- 
terial is weak to moderately correlated |6|, since specific 
heat enhancements are moderate and no convincing satel- 
lite or Kondo photocmission peaks are observed, which is 
consistent with the success of band-structure in predict- 
ing materials properties 0, @|. At the same time, when 
the uranium atoms are pushed apart by other elements, 
they form many heavy fermion and other strongly cor- 
related uranium compounds 9] . In this regard, uranium 
is an inviting target to study, since it should have in- 
teresting correlation effects beyond conventional metals 
like copper or aluminum, and yet these should be weak 
enough to have some hope of accurately calculating them. 
It is thus an important testing ground for correlation the- 
ory and how many-body effects correct conventional LDA 
band structures. 

The most widely used electronic-structure method, the 
local density approximation (LDA), has been an im- 
mensely successful tool that reasonably predicts ground- 
state properties of weakly correlated systems. The LDA 
is much less successful at predicting optical properties of 
such systems, and its failures become more serious as cor- 
relations become stronger. Recent photoemission spec- 
troscopy on high quality uranium singe crystals has re- 
vealed additional information about the electronic struc- 



ture of this material [lCj, lll|. Comparison with LDA 



calculated electronic bands shows some disagreement be- 
tween experiment and theory. Because of the poor treat- 
ment of electron correlations by LDA it is difficult to 
conclude whether the observed discrepancies between the 
predicted band structure and photoemission data are due 
to electronic correlations, even though perhaps weak, or 
to other effects such as surface states. For the same rea- 
son it is not clear how much of the mass enhancement 
observed in the specific heat [3, 12 1 can be attributed to 
electron correlations and how much to electron-phonon 
coupling. To date all first-principles theoretical treat- 
ments of the uranium electronic structure have been 
based on LDA or the generalized-gradient approximation, 
GGA, extension to LDA. Therefore, it is important to ex- 
plore the electronic structure of uranium with methods 
that treat more accurately the electron-electron interac- 
tions and to understand how they affect the electronic 
properties of this material. 

In Pu, it is now standard to use dynamical mean-field 
theory (DMFT) to treat the strong correlations that go 
well beyond conventional band-structure However, 
this has the unsatisfactory aspect that a model Hamilto- 
nian is grafted onto a band-structure approach in an ad 
hoc manner. Because of the much weaker correlations in 
uranium, it is possible in this material to instead use a 
more approximate treatment of correlation effects that is 
completely first principles and yet goes significantly be- 
yond conventional band theory. Thus, in this Letter we 
apply for first time for any actinide a rigorous ab initio 
self-consistent many-body theoretical approach, the GW 
approximation, and show how electronic correlation ef- 
fects modify the electronic structure of uranium that is 
predicted by LDA band theory. 

In this work we use the recently developed QSGW ver- 
sion [li , HI, IH of the GW method, which itself can be 



2 



viewed as the first term in the expansion of the non-local 
energy dependent self-energy S(r, r, u>) in the screened 
Coulomb interaction W. From a more physical point 
of view the GW approximation can be interpreted as a 
dynamically screened Hartree-Fock approximation plus 
a Coulomb hole contribution [T^. Therefore, GW is a 
well defined perturbation theory. In its usual implemen- 
tion, sometimes called the "one-shot" approximation, it 
depends on an one-electron Green's functions based on 
LDA eigenvalues and eigenfunctions, and the results can 
depend on this choice. We have demonstrated [l5| that as 
correlations become stronger serious practical and formal 
problems can arise in this approximation. In Ref. [3] a 
formal description is provided on how QSGW is a rig- 
orous way to surmount this difficulty, based on using a 
self-consistent one-electron Green's function that is de- 
rived from the self-energy (the quasi-particle eigenvalues 
and eigenfunctions) instead of an LDA starting point. In 
the literature, it is has been demonstrated that QSGW 
version of GW theory reliably describes a wide range of 
spd systems 3, 18, HI and rare-earths [5(|. It should be 



noted that the energy eigenvalues of the QSGW method 
are the same as the quasiparticle spectra of the GW 
method (i.e., the peaks in the self-energy). This captures 
the many-body shifts in the quasiparticle energies. How- 
ever, when presenting the quasiparticle DOS this ignores 
the smearing by the imaginary part of the self-energy 
of the spectra due to quasiparticle lifetime effects, which 
should increase as one moves farther away from the Fermi 
energy. 

The QSGW method is currently implemented using a 
generalization of the Full Potential Linear Muffin Tin Or- 
bital (FP-LMTO) method [2l|, so we make no approxi- 
mations for the shape of crystal potential. The smoothed 
LMTO basis [IB} includes orbitals with I < l max — 6; both 
7p and 6p as well as both 5/ and 6/ are included in the 
basis. 6/ are added in the form of local orbitals [HI], that 
is an orbital strictly confined to the augmentation sphere, 
and has no envelope function at all. 7p are added as a 
kind of extended local orbitals the 'head' of which is eval- 
uated at an energy far above Fermi level [l5j], instead of 
making the orbital vanish at the augmentation radius a 
smooth Hankel 'tail' is attached to the orbital. A partic- 
ularly important point is that core states are treated at 
the exchange-only level. We have demonstrated in some 
detail ID, 22| that approximating the core by the LDA 
potential, i.e. computing E from the valence electrons 
only, sometimes leads to significant errors. Since QSGW 
gives the self-consistent solution at the scalar relativistic 
level, we add the spin-orbit operator Hso — L ■ S/2c 2 as 
a perturbation (it is not included in the self-consistency 
cycle). For our calculations we use the equilibrium crys- 
tal structure of a-U, the orthorhombic Cmcm, with the 
Uranium atoms located at the 4c positions:(0,y, j) and 
(0,-y,§) plus C centering; we use the experimental lat- 
tice parameters a = 2.858A, b = 5.876A, c = 4.955A, 
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FIG. 1: The energy bands (or quasi-particle energies) along 
two symmetry directions (left panel) and the total DOS (right 
panel); the QSGW results are represented by solid red lines 
and the LDA band-structure results by dashed blue lines. 



and y=0.105. 

In Fig. 1 we compare the calculated QSGW one- 
particle electronic structure of a U with the LDA band- 
structure results. In both cases, the narrow bands located 
approximately between -1 and +3 eV are mainly due to 
the uranium 5/ orbitals; the lowest dispersive bands seen 
on this plot have s character; and the bands above 3 eV 
are strongly hybridized. Both methods are roughly in 
agreement for the total density of states (DOS). How- 
ever, the large DOS peak that is a little above the Fermi 
energy, Ep, is narrower in width and larger in magnitude 
for the QSGW calculation; also, the quasiparticle ener- 
gies only agree well with the LDA band-structure results 
in the vicinity of Ep. As we move to higher or lower ener- 
gies (away from Ep) the difference between QSGW and 
LDA quasiparticle energies gradually increase. Among 
the occupied states, the metallic s-d bands at the lowest 
energies experience a significant downward shift relative 
to the / bands when compared to the LDA results (note 
that the main part of the p states are believed to lie 
above Ep, since they are repelled by the 6p semi-core 
states that fall below and are well separated from the 
conduction band, and hence only p hybridization tails 
appear in the occupied conduction-band region). For ex- 
ample, at the r point the shift is about 1 eV downwards, 
but more generally, however, the energy shift is some- 
what k-dependent. The partial DOS presented in Fig. [3] 
shows that after integration over all k there is a down- 
ward energy shift of about 0.5 eV for the occupied s-d 
bands. 

In Fig.[^b) we present the partial DOS for the / bands. 
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FIG. 2: (a) The d partial DOS (b) The / partial DOS (c) 
The s partial DOS. In all cases the QSGW results are the 
solid red lines, and the LDA results are the dashed blue lines. 



One of the important effects of electronic correlation is to 
narrow the width of a band. This shows up as a narrower 
/-band width in the QSGW calculation. In addition, 
since the area under the curve is proportional to the num- 
ber of / states, which remains constant, the amplitude 
of the quasiparticle peak is also higher. The narrowing 
of the / band together with its energy shift results in a 
slight change of the electron occupation. Comparison of 
the partial DOS shows that the / band shifts up, rela- 
tively to s and d bands, toward unoccupied states. The 
/ occupation is equal to 3.19 in the QSGW and 3.57 in 
the LDA calculation, hence in QSGW about 0.4 / elec- 
tron is lost to the s-d interstitial charge. The overall / 
occupation in uranium is relatively low so that the posi- 
tion of the Fermi level remains in the lower part of the / 
peak, where the difference between QSGW and LDA cal- 
culation is negligible. For this reason, even though the 
5/ electron states in uranium appear to be correlated, 
the physical properties that are related only to occupied 
electron states should be predicted well by the LDA ap- 
proximation. It is mainly in the excited-state spectra of 
the / states above the Fermi energy where the correla- 
tion effects are strongly apparent. Consistent with the 
weak to moderate strength of the correlations, we find 
that the first iteration of the QSGW method, which is 
sometimes called the "one-shot GW" is very similar to 
the fully self-consistent QSGW results. 

The electron DOS at the Fermi level in QSGW is 
3.35 states/(eV atom) while in the LDA result is 2.75 
states/(eV atom). We have found that the DOS has a 
lot of variation around the Fermi level and requires a very 
large number of k points to converge (we used 82 3 ). This 
explains deviations of our DOS with earlier results (e.g., 
Ref. [llj and references therein) Our LDA LMTO results 
using a Barth-Hedin exchange-correlation potential were 
found to be in excellent agreement with GGA (gradient- 
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FIG. 3: (a) Comparison of photemission data with the calcu- 
lated QSGW DOS (b) Comparison of photemission data with 
the calculated LDA DOS. In both cases, the DOS has a Gaus- 
sian broadening of 0.05 eV 



corrected) LAPW calculations using a PBE exchange- 
correlation potential that we calculated with the Wien2k 
code [23j | : e.g., we found insignificant differences in the 
total LDA and GGA DOS. 



It was shown by Luttinger [24( that the linear coeffi- 
cient of specific heat for a system of interacting electrons 
is given by the quasiparticle DOS, i.e., 7 ~ ^2 k S(EF ~ 
E(k)). This can be compared to the one-electron coef- 
ficient, which is given by 7 ~ J2k^(-^F ~ e (k))- Here, 
Ep and E(k) are the Fermi and quasiparticle energies 
of the interacting electron gas (in the QSGW approxi- 
mation), while E F and e(k) are the Fermi energy and 
band-structure eigenvalues. Hence only the quasiparti- 
cle shifts are needed for calculating the specific heat, 
which are included by construction in the QSGW en- 
ergy eigenvalues and DOS. Therefore, the specific heat 
of the interacting QSGW electron gas is proportional to 
the QSGW DOS at the Fermi level N{E F ), and there 
is no need to include an additional renormalization fac- 
tor (1 — dYi/du)) factor, which in model calculations con- 
verts the band-structure DOS to the quasiparticle DOS. 
We find that the linear coefficient of specific heat in 
QSGW is 7= 7.89 mJ mob 1 KC 2 while in LDA is 7= 
6.48 mJ mol -1 K -2 , giving a QSGW enhancement factor 

Of NQ SGW (E F )/N LDA {E F ) = jQSGW/^LDA = 122 

A recently measured value (t2| . r y exp =9.l5 mJ mol -1 
K~ 2 , is larger than the QSGW result by a factor of 
lexp/lQSGW— 1-16 and larger than the LDA result by a 
factor of 7 e a:p/7LDA— 1-41. The relatively small remain- 
ing enhancement that is not accounted for by QSGW 
must almost certainly be the electron-phonon enhance- 
ments that are present in all metals and that are typically 
at least this big. 
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FIG. 4: The Fermi surface cross section cut by the [100]- [010] 
plane at h=0.288; QSGW - solid lines, LDA - dotted lines. 



The bremsstrahlung isochromat spectroscopy (BIS) of 
a uranium exhibits a double peak structure within the 
interval of 1 to 2 eV above Ep [25]; as one can see in 
Fig. [3J this feature is reproduced by both calculations. 
However, further comparison of the two methods with 
photoemission data is difficult because the photoemis- 
sion peaks are too broad comparing to the band width 
narrowing and quasiparticle shifts observed with the help 
of QSGW method. Recent angle resolved photoemission 
(ARPES) [ll[ and ultraviolet photoemission spectroscopy 
(UPS) [10|| were found to be in good agreement with GGA 
band structure calculations. However, a low energy UPS 
peak and several APRES local maxima were not pre- 
dicted by GGA band structure. Since the QSGW energy 
bands along r — ► Z direction in the energy window -2 
to eV are very similar to those obtained with LDA the 
agreement with the valence band (UPS) is of the same 
level with LDA. The unexplained APRES local maxima 
are located in the vicinity of T point with energies around 
-5 eV and -2 eV 0. Even though the QSGW bands 
along r — > S direction agree less with LDA than they do 
along r — > Z direction, the differences are not significant 
enough to explain the aforementioned ARPES data. We 
suspect that the presence of uranium surface states as 
was suggested in Ref. [ll[ or final state effects are a more 
likely explanation for these features. 

The [100]- [010] cross section of the Fermi surface is 
presented in Fig. 3. It shows that there are only slight 
changes in QSGW compared to the LDA calculations. 
This is representative of several cross sections that were 
plotted, showing that overall the electron correlations 
have no significant effect on the shape and size of the 
Fermi surface. 

In conclusion, we have used the QSGW method to 
show that moderate /-electron correlation effects arc 
present in cv-U, and that it is because of the low oc- 
cupation of / electrons that these effects don't show up 
more strongly in Fermi surface and other ground-state 



properties of this material. Most of the correlation ef- 
fects only appear in the excited-state spectra in the un- 
occupied / states. As is commonly suspected, LDA or 
GGA band-structure methods somewhat misplace nar- 
row bands, such as d and / bands, with respect to the 
remaining metallic bands. For uranium the error is about 
0.5 eV. 
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